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Abstract 

A new nonparametric estimator of the local Hurst function of a multifractional Gaussian process based on the 
increment ratio (IR) statistic is defined. In a general frame, the point-wise and uniform weak and strong consistency 
and a multidimensional central limit theorem for this estimator are established. Similar results are obtained for a 
refinement of the generalized quadratic variations (QV) estimator. The example of the multifractional Brownian 
motion is studied in detail. A simulation study is included showing that the IR-estimator is more accurate than the 
QV-estimator. 

Keywords Nonparametric estimators Hurst function tangent process multifractional Brownian motion Gaussian 
process central limit theorem. MSC Primary: 62G05 Secondary 62G20 60F05 60G22 

1 Introduction 

The present paper discusses nonparametric estimation of the local Hurst function H(t),t € (0, 1) of a multifractional 
Gaussian process Z = (Z{t)) te ^ Q ^ observed at discrete times k/n, k = 1, • • ■ ,n — 1. The process Z is characterized by 
the fact that its infinitesimal increments at each point t € (0, 1) resemble a multiple of a fractional Brownian motion 
J3y( t ) with Hurst parameter H (t) <E (0, 1): 

{\- H ^{Z{t + \r)-Z{t))) T ^ {c(t)B m (r)) T > , (1.1) 

in the sense of the weak convergence of finite-dimensional distributions, where c(t) > 0,t 6 (0,1) is a deterministic 
function. Recall that a fractional Brownian motion (FBM) with Hurst parameter H 6 (0, 1) is a Gaussian process 
(Bh(t)) t>0 having stationary increments and satisfying EBh{t) — and EB|(t) = t 2H ,t > 0. Property (jl.lj) is 
called asymptotic self-similarity and the limiting process a tangent process at t. For a Gaussian zero-mean process Z, 
relation (jl.lj) is equivalent to 

E(Z(t + r) -Z(t)f ~ c 2 {t)T 2H ^\ r^0. (1.2) 

A typical example of such Gaussian process Z is a multifractional Brownian motion (MBM) defined through the 
following harmonizable representation 

X(t)= B H{t) (t) = K{H{t)) J ^~ + \ /2 W{dx), with K{H):={HT{2H)^H)/k) 1, \ (1.3) 

where W{dx) is a complex- valued Gaussian noise satisfying E[J M g(x)W(dx)J M h(x)W(dx)] = J M g(x)h(x)dx 7 h,g e 
L 2 (M.) (sec more details in section 7.2.2 of [H]) A first version of MBM was first introduced by Peltier and Levy- 
Vehel 20J using a slightly different time-domain representation, while the MBM in (| 1 . 3[) was defined, studied and 
generalized in several papers. In particularly, Stoev and Taqqu |22j defined a generalization Yf a + tCt -\ of p.3j) depending 
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on parameters (a + ,a~) £ M 2 \{(0, 0)} such that ^(1.1) agrees with X of (|1.3[) up to a multiplicative deterministic 
function, see Section 0] below for details. In the present paper, the functions H(-) and c(-) are assumed to belong to a 
Holder class C(0, 1), 7/ > 0. The smoothness parameter r\ plays an important role in the asymptotic results discussed 
below. In particular, for < r\ < 1 the local Hurst function of MBM can be very irregular making the problem of its 
estimation more difficult. 

Another important parameter for estimation of H{t) is the bandwidth parameter a £ (0,1). Following a general 
approach in non-parametric estimation, discrete observations Z(k/n), 1 < k < n — 1 are localized in a neighborhood 
V nt a.(t) = {k/n : \k/n — 1\ < n~ a } of a given point t £ (0, 1) containing at most [2n 1- "] + 1 points, and an estimate of 
H(t) from the localized sample Z(k/n), k £ V n , a (t) is constructed. Choosing a smaller a means that the estimate of 
H(t) uses more observations of Z and therefore the variance of it can be expected to be smaller. On the other hand, 
increasing the "width" of V n ,a{t) (by taking a smaller a) usually increases the bias of the estimator particularly when r\ 
is small, since the function H (•) is more likely to change in a larger interval V„ jCt (i). Hence the accuracy of estimation 
of H (t) depends on a and 77 in a crucial way. 

The present paper discusses two discusses two estimators of H(t). The first estimator Hn^} (t) is a localized version of 
the increment ratio (IR) estimator introduced in [5] . It is written as a localized sum of ratios of second-order variations 
V£ Z(k/n) — Z((k + 2)/n) — 2Z((k + l)/n) + Z(k/n), k £ V n ,a(t) and is related to zero intersections of the variations 
sequence, see (|2.8|) . also ([5,, Remark 1.1). We also discuss a new estimator Hn,a 2 \t) which is a pseudo-generalized 
least squares (PGLS) version of (t) constructed from several dilatations of the second-order variations. At the 

same time, we study the generalized quadratic variations (QV) estimator H^a\t) of [15] and its refinement (a PGLS 
version H$P\t). The main reasons for including the QV estimators in this paper is the fact that the study and 
the asymptotic results for IR and QV estimators are very similar (both estimators are expressed in terms of sums of 
nonlinear functions of Gaussian vectors of the Hermite rank 2), and our desire to compare finite-sample performance 
of these two classes of estimators. 

Let us briefly describe the main results of this paper. Section [2] introduces the two main Assumptions (A) K and 
(B) a on a Gaussian process Z . The first assumption specifies the rate of convergence in a similar asymptotic relation 
to (|1.2[) and controls the bias of our estimators. The second assumption concerns the decay rate of the covariance of 
generalized variations of Z and thus helps to control the variance of the estimators. Under these assumptions, Theorem 
Q] (i)-(iv) obtains the weak and strong consistency of the estimators H^a\t) and H^a \t) at a given point t £ (0, 1) 
and some uniform weak and strong consistency rates of these estimators over arbitrary interval t £ (e, 1 — e), e > under 
increasingly restrictive conditions on a, rj 1 k and other parameters. A multidimensional central limit theorem (CLT) 
satisfied by the above estimators is proved in Theorem [2] Corollaries Q] and [2] extend Theorems [T] and [2] for the PGLS 
versions of the IR and QV estimators which are studied in Section [3] Proposition [3] provides limit theorems to the IR 
and QV estimators for the MBM (including its generalization in [22]). Section [5] is devoted to simulations while Section 
[S]and Section [7] contain the proofs. Let us note that the simulations in Section [5] performed on the MBM in (jl.3[) for 
some typical Hurst functions and three different values of the smoothness parameter 77 show a very good performance 
of the estimator Hn,a 2 \t). The last estimator provides in almost all cases a smaller MISE criterion than (t). 
We also note that the numerical results for both classes of estimators are clearly (in the case of IR) or slightly (in the 
case of QV) in favor of their PGLS versions. 

Nonparametric estimation of the local Hurst function was the subject of several papers (pQ, [3], [8], [9], [10], [TT] . 
[To] . |18j . [19]). However most of these results are less comprehensive and sharp when applied to the MBM while some 
contain erroneous statements. See Remarks [21 El and [51 below, for more detailed discussion and some controversies. 
The MBM is probably the most important example of a multifractional Gaussian process at the present time. It is 
considered as a benchmark for the estimation problem of H(t) in the present paper and is discussed in detail in Section 
[H We believe that the conditions of Proposition [3] for the IR and QV estimators of the Hurst function of the MBM 
are optimal or close to optimal since they are derived from the asymptotic expansion of the covariance of generalized 
variations of the MBM in Lemma [U See also Remark [2] 
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Notation. In what follows, C stands for a constant whose precise value is unimportant and which may change from 
line to line. Also, we write — ^> , , — ^> for convergence in probability, almost sure (a.s.) convergence and 

n— ►oo n— ¥00 n— >oo 

the (weak) convergence of probability distributions, respectively. [77] denotes the integer part of a real number 77. 
C v (a, b) is the space of all [77] > times continuously differentiable functions / : (a, b) 1— > R on interval (a, b) C R such 
that /CM), the [77] th derivative of /, is an (77 — [77] )-H61derian function on (a, b), i.e. there exists K > such that 
\f(li])( x ) - /(W)(y)| < if I x -y l^-W for all (a, 6). 



2 General results 

2.1 Definitions of estimators 

Consider a filter a := (ao, • • • , a g ) € R 9+1 such that there exists m £ N* satisfying 

^~Va^ = for p = (),•■• ,(m-l) and ^^ m a £ ^0. (2.1) 

f=o 1=0 

Denote A{m, q) the above class of filters. For a £ „4(777, g) and 71 £ N* define the corresponding generalized variations 
of Z by 

V*Z(t) := J2 a i z (t + *M, V a Z(t) ee V?Z{t). (2.2) 

e=o 

The main examples of generalized variations are the usual simple variations corresponding to a = (1,-1) £ -4(1, 1), m = 
q = 1, and the second order variations corresponding to a — a* := (1,-2,1) 6 .4.(2,2), m = q = 2. For a filter 
a = (ao, • • • , a ? ) £ A(m, q) and j€N define its jth dilatation = (a^ , • • • , a^) € A(m,jq) by 







W) — „. o,^ — 



a, and a^ ; := if fc <^ jN. 



For a € (0, 1), define a neighborhood of t £ (0, 1) and its cardinal by: 

V n ,a(t) := {fc € {1,2,- •• ,n-g- 1}, |- - 1| < n""} and v n , a (*):=# V„,o(t). 

Note that for any i G (0, 1) and n large enough, u„. Q (t) = [2n 1 ~ Q ] or [2n 1_Q! ] + 1 depending on the parity of [2n 1_Q ] = 
the integer part of 2n 1 ~ a . For H £ (0, 1), define 

A 2 (H) := E^(vfB H (0),VfB H (-)), (2.3) 

71 

where 

*(* (1) .* (a> ) - pSjrA' (2.4) 

is a function taking values in [0, 1], more explicitly, 



MH) = - arccos(- P2 (ff)) + -J - + ^ log ( 2 ), (2.5) 

p 2 (ii) := Cor(Vf B H (0),V n a ' B H (-)) = g _ 22g+1 . (2.6) 

The function A2 does not depend on n by self-similarity of FBM and is monotone increasing on (0, 1), see [5]. The 
estimators H^a(t) and Hn$(t) are defined for t £ (0, 1) from observed sample (Z{l/n), Z(2/n), • • • , Z((n — l)/n)) 
as follows: 

: =^M^ E <"> 

' v ; fceVn, a (t) — 



\« n a [t) * — ' 77 Tl J 



(2. 
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where A := (logi — ~ X)j=i ^°SJ')i<i< ^ M p is a row vector, A T its transposed vector (vector-column), and A 2 1 is the 
inverse function of A 2 . Thus, H^a^ if) is defined from a log-regression of generalized quadratic variations, see also [15] . 
whereas Hn.a^ (t) is defined from a sample mean of ratios of the second-order variations in a neighborhood of a given 
point t £ (0, 1). A discrete time version of the last estimator was first studied in [24] . 

2.2 Assumptions 

First, recall (see e.g. [24]) that for a filter a = (do, • • • , a g ) £ .4(m, g) with m > 1 and for any fixed H £ (0, 1) and 
n £ N, (V^ Bff(j/n))j eZ is a stationary Gaussian process, (V^-Bff(^)) . gz ^= d ' n~ H (V a BhU)) j £Z by self-similarity, 
furthermore, 



Var(y n aU J B ff (i))=-2 F C(iT ! a), C(tf, a) := Var(V Q B ff (0)) = -- ^ a fc a fe , |fc - fc| 2H , i€N', 

k,k'=0 

Cov(v:B H (^V a B H (^))=-^jrj2^ki-k + j- 1 \ 2H 

i=0 j=0 

/ (V 9 a i m ) 2 2 ^ \ Ifr' - k\ 2H ~ 2m 
~ W II ^ - £ )) . - fc l - ». (2.9) 

Note C(H,a) > for any H £ (0,1), a £ A(m,q) which fact follows from the harmonizable representation of FBM 
analogous to (|1.3[) yielding C(H,a) = K 2 (H) J R \ J2l=o a k{ elhx ~ l)\ 2 \x\~ 2H ~ 1 dx > since the integrand does not 
identically vanish on R. In the particular case a = a* = (1, —2, 1) (m = 2), then 

Vex(V a * B H (j)) = 4-4 H , 
Cor(V a *B H (j),V a *B H (j + l)) = P2 (H), 

where p2(H) is defined in (|2.6I) . 

For a given m > 1, a general Gaussian process Z = (Z(t))te(o,i) and real numbers n > and a £ (0, 1), we introduce 
the following assumptions: 

(A)« There exist continuous functions < H(t) < 1 and c(t) > for t £ (0,1) satisfying H(-) £ C(0, 1) and 
c(-) £ C(0, 1), r\ > such that for any < e < 1/2, j > 1 and any filter a £ A{m, q), q>l, 

V&v(V^Z(k/n)) fe , . 

V " W " -c(-) — > 0. (2.10) 



max n 

H<K[(l-e)n] 



Var (V-B H(k/n) {0)) V 



(B) a There exist C > 0, 7 > 1/2 and < 9 < 7/2 such that for any n £ W, < e < 1/2, [ne] < k < k' < [(1 - e)n] 
and any filter a £ A(m, q), q > 1, 

Cor(F n Z(-),F n °Z(-))| < Cn {1 - a)6 {\k' -k\ An 1 -")" 7 . 



Remark 1 If condition (A)o is satisfied for m = 1 then Z admits a FBM as a tangent process at each point t £ (0,1), 
in the sense that for any k £ Z, 

n^)(z(M±*±i)_ Z( M±*)) ^ (^ ff(t)(i) ) . (2.11) 

Assumption (A) K with k > allows to control the uniform convergence rate of generalized variations of Z to the 
corresponding variations of the tangent process. Condition (B) Q provides a bound for the correlation of the process 
(V£Z(k/n)) r„ e ] <fc fc/<r(i_ e ) n i ( tne choice of this bound can be deduced from the expansion (|2.9[l ). Under such conditions, 
the strong consistency and a CLT for the estimators Hnfa\Hn$ respectively defined in (|2.7|) . (|2.8|) can be established 
(see Sec. 2.3 below). 
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The following statement deduced from Assumption (A) K will be often used in the remaining sections. 
Property 1 Assumption (A) K implies that for any a € A(m,q) and any s > 0, j S N 



max n 

[ne]<k<[(l-e)n] 



Coi{v:Z(^),V:Z(^))-Cor(V a B H{i) (0),V a B siA) (u)) -> 0. (2.12) 

77. 77. n ' TL ' iri rv, 



Proof. We have 



Hv:z(t)v:z( k -±i)] = l{^(v:z^)+v^z(^))-v a v(v n ^z(l)-v:z( k -±i))} 

= i{Var(^ + ^))-Var(^-Z(^))}, 
where generalized variations Z(^) = Ya=o a i = -^(^r) correspond to filters a ± = (a*, • • • , a^ +J ) S with 

at := Oil(0 < i < q) + ai-jl(j <i<j + q), a- := Oil(0 < i < g) - < i < j + g). (2.13) 

Note that a € _A(m, g) implies a € A(m,j + g) for any j = 0, 1, ■ ■ • . From (A) K with a — a we have that uniformly 
in [ne] < k < [(1 - e)n], V n Q± Z(|) - c(|)VarF n 0± S Jf(fc/n) (0)(l + o(n-«)), implying 

Cov(v; a Z(^),^Z(^±2)) = ic(J){Var(F„ a+ J B H ( fe /n)(0))(l + O (ri- K )) - Var(y-~B H(k/n) (0))(l + O (n-"))} 

= ± c (£){Vax(V? + B ff(fc/n) ^ 

= C (-)Cov(V„ Q i^ (fc/n) (0)^^ (2.14) 
n n 

We also have 

K a Z(-) = c (-)Var(y n a B ff(fc/n) (0))(l + O (n-«)), ^Z(^) = c(-)Var (^5 ff{k/n) (0))(l + o( n -*)), (2.15) 

'Ai TXi 7~b 

where the last relation follows by writing Var (V£Z = Var (V„ Q+ Z(£))-Var (F n a Z(|))-2Cov (V„ a Z(£), V£Z{^-)) 
and using the facts established above. Then, (|2.12p is deduced from (|2.14[) and (|2.15p . □ 

2.3 Asymptotic results 

Denote 

2(l-a)(l-i)-a 2(l-a)(l-f)-a-l 

/x:= ^ , /i i: = ^ . (2.16) 

Note /i > (respectively, pb\ > 0) is equivalent to a < |^E§g (respectively, a < X~ 2 2 g )- 

Theorem 1 Lei Z = (Z(i)) te ( ,i) ^ e a zero-mean Gaussian process. 

(i) If Z satisfies condition (A)o and (B) Q wif/i m = 2, a £ (0, 1) then for any t G (0, 1) 

(w^ // Z satisfies conditions (A)o and (B) a m = 2, < a < 2 |~^ £/zen /or any t G (0, 1) 

fmj //Z satisfies conditions (A) K and (B) a m = 2, k > fi and 3j _2o+4y( 7] /\2) — a ^ 37-20 ^ eri ■ / " or e > 

sup |5^(t)-i?(t)| = O p (n-"). 

e<t<l-e 

fivj // Z satisfies conditions (A) K and (B) Q with m = 2, re > and 3 7 _ 2 e'+4^( T) A2) — a ^ §7^9 i/ien /or any 
e > 0, (5 > 

sup |il^(t) -/?"(*) I = 0(n-^- 5 )) a.s. 

e<i<l-e 
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Theorem 2 Let Z = (-Z(i))te(o,i) be a zero-mean Gaussian process satisfying (A) K and (B) Q , with m = 2, a > 
i+2(i)A2) ' K — an d ^ — 0- TTien /or any n£N* and any t 1; • • • ,t u £ (0, 1), t, 7^ ij (i 7^ j) 



<i< t 



(w^(*i)) a 



<i<i 



where W^ IR \ti), i = 1, • • • , it are independent centered Gaussian r.v. 's such as, with A2(-) defined in i2.5\) . 
E[W^ R Ht t )r := [^(AaJ-^Aa^t,-)))] 2 ^^)) 



(2.17) 



(2.18) 



a 2 (tf):=^Cov ^(r'BH^r'B^l^fr'^W.r'Bfftfc + l))) /orffe (0,1). (2.19) 



The following Theorem 3 extends Theorems Q] and [2] to the estimator H^f ^ . As noted in the Introduction, both 
estimators H™{t) and are written in terms of sums of nonlinear functions of Gaussian vectors of the Hermite 

rank 2 and their asymptotic analysis is very similar. The proof of Theorem 3 can be deduced mutatis mutandis from 
the proofs of Theorems Q] and [2] and is omitted. 

Theorem 3 Let Z be a zero-mean Gaussian process and a G A(m, q), m > 1 a filter. 

(a) The statements of Theorem]]] (i)-(iv) hold for H^a\t) instead of Hn,a\t) under respective assumptions (i)-(iv) 
with m = 2 replaced by m > 1. 

(b ) Let Z satisfy the assumptions of Theorem [D with m = 2 replaced by m > 1 . Then for any u £ N* and any 
tl,»- ,t u G (0,1), U^tj {i^j) 



yfi^mWM *fc))i<,<„ ^ {w {QV) tU)) 1<i<u , 

— — 71— »0O — — 

where W^ v '(ti), i — 1, • • • , u are independent centered Gaussian r.v. 's such as, with A defined in f£. 



(2.20) 



V[W iQV) (U)} 



A T T(H(ti))A 



and r 



2H_ 



v ' 12 jez 2^*;i,fe 2 =o °fei a fc 2 l Kl K2 I 

3 Pseudo-generalized least squares versions of the IR and QV estimators 

Asymptotic and finite-sample performance of the estimators Hn$ (t) and H^a^ (t) can be improved by using their 



l<ii ,i2<p 



,(2.21) 



pseudo-generalized least squares (PGLS) versions Hn.a 2 \t) and Hi^J z, (t) as defined below. 

Definition of H [ r ™\t). Let a^*~> denote the ith dilatation of the filter a* = (1, -2, 1), i e N*, so that V^"' Z(£/n) 
Z{{1 + 2i)/n) - 2Z((£ + i)/n) + Z(£/n). The corresponding IR estimator of H(t) is defined as 

where the function A^ : (0, 1) — > (0, 1) is defined as the expectation 



HQV2), 



A^(H):=E^(vf' : 'B H (0),vf ' B„(-)) = - arccos(-^ ) (H)) + - 



n 7r 



loe 



and ^ m :=C, K "" Bl (»,,,r*(i)) - + H" - l«- H" + 4|i - 



Under the assumptions of Theorem [21 these estimators satisfy the p— dimensional CLT: 



n (l-a)/2 (HW>*(t) - H(t) 



l<i<p n— too 



M p (0 , EW(fl-(t)) 



(3.1) 
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with the limit covariance matrix EW(iJ(t)), where for H g (0, 1) 

^ {P){H) ■= (II^^ ^^^ (3 - 2) 

' hi \\V a ™B B $)\ + \V*™B H {V)y \V°™B H {k)\ + \V« w B H {k + l)\) [ ' 

(note that V 0<i *' B H (k) = B H (2j + k) - 2B H (j + k) + B H (k) for j = 1, • • • ,p and fc € Z). Now define 

E«(t) := (3.4) 

which is a consistent estimator of E( p )(£f(i)) since Hn^ (t) is a consistent estimator of Hit). Then, with the vector- 
column l p = (1, 1, • • • , 1)t g MP . the PGLSE (t) of H(t) is defined by 

H { Z 2) {t) := (IT (EW(t))- 1 lp)" 1 IJ (EW(i))- 1 . (3.5) 

Definition of Hj^ 2 \i). Analogously to ([3. 4ft , define 

f(t) := r(H<M{t)), (3.6) 

which is a consistent estimator of T(H(t)) since is a consistent estimator of H(t). The one-dimensional CLT 

satisfied by H^a\t) at a given point t g (0, 1) (see (|2.20p ) is a consequence of the following p— dimensional CLT: 

V2^(f$(t) - 2ff(t) (log(t/n)) - C(i)lp) A 7V P (0, r(H(t))) 
with := log ( - c -f n,*'=o " fc l 2H(t) ) and 

see [B] for details. Then 2H^a\t) is the slope of the ordinary least squares regression of the vector T^) a (t) onto 
the vector ( log(i/n)) 1<i<p . The PGLS estimators H^ 2 \t) and ci Q Q V2) (t) of £T(t) and C(t) are obtained from the 
corresponding generalized least squares regression with the covariance matrix T(H(t)) replaced by its estimate T(t) 
(EH): 

( 2 |f^w) := (^ )F W" 1 (^ ) ) T )~ 1 ^ )F W" 1 ^w) T ' (3 - 7) 

with 4 ri := ^ 1 ° g( i 1/n) ^/n) ... log(p/n)Y (3.8) 

Corollary 1 (jj Assume that for each < H < 1 the matrix T,^(H) of \3.2\) is non- degenerate and continuously 
depends on H. Then all statements of Theorem]]] (i)-(iv) hold for the estimator ifj^a instead of Hnfa 
(ii) Assume that for each < H < 1 the matrix T(H) in h2.21}) is non- degenerate and continuously depends on H . 
Then the statement of Theorem 3 (a) holds for the estimator H^a^ instead of H^f^ ■ 

Corollary 2 Let Z satisfy the assumptions of Theorem^ and the matrices E^(iJ), T(H), < H < 1 satisfy the 
conditions of Corollary^ (i)-(ii). Then for any hgN* and any ti, ■ ■ ■ ,t u € (0, 1), ij ^ tj (i ^ j), 

^^(HLT^)-H(U)) 1<l<u A (W^(U)) 1<i<u (3.9) 

— — n—tao — — 

and 

In Iff.ffp . W^ V2 ^ (ti), 1 < i < u are independent centered Gaussian r.v.'s with respective variances T^[W^ IR2 \ti)] 2 := 
(lj (E^(fl r (tj))) _1 l p ) . In \3.10}) , W^ V2 \ti), 1 < i < u are independent centered Gaussian random vectors with 
respective 2x2 covariance matrices E^W 72 )^)^^ 2 )^)^ := (z[ p) r^^))" 1 (^ p) )t)~\ 
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4 The case of mult ifr actional Brownian motion 



This section details Assumptions (A) K and (B) Q for the MBM (|1.3|) and its generalization due to Stoev and Taqqu 
[22] . As a consequence, the asymptotic behavior of the IR and QV estimators for the above process is established under 
explicit conditions on the parameters a, ?y and the function H(-). 

The above mentioned generalization of the MBM in [22] is given by stochastic integral representation 

T(H(t) + -) f e itx - 1 

*(«,+,„-)(*) := 7= 2 / , , gm+ i U ia+ia - ) (H(t),x)W(dx), feR, (4.1) 

V 27T JR |X| w " t *2 

where (a + , a - ) € IR 2 \ {(0, 0)} are parameters, W^da;) is the same white noise as in (|1.3[) . and 

:= a + e- isi S n ^( H +^f + a^e'^^^+^f , H e (0,1), z e R. (4.2) 

The above definition implies EF ( 2 + .(i) = EF ( 2 + Q - ) (l)|i| 2ff(t) and a covariance structure of F( a + a -) similar to that of 
the MBM in (|1.3j) (see [22], Theorem 4.1). In particular, Yr\ $\ gives a harmonizable representation of the time-domain 
MBM of [20]. The so-called well-balanced case case a + = a~ = a ^ 0, or 

w) = ^cosa^+i)^^^)/^^^) 



is more delicate, see ([22], sec. 5), since the function H i-> cos ({H + |)S ) vanishes and changes its sign a.t H = 1/2. In 
particular, F( 0i0 )(t)/(EF ( 2 J . q) (1)) 1/2 = -sign(F(t) - where A is defined at fOjl . Since X has a.s. continuous 

trajectories under the condition (|4.5|) on H(-), see p], this implies that the normalized process F( Qia )(t)/(EF^ 2 ^(l)) 1 / 2 
is generally a.s. discontinuous at each t ^ with H(t) = 1/2. 

In order to avoid the above complication, we introduce the following generalization AV a + )Q -\ of the MBM in (|1.3|) : 

Y M - ! a+=a-^0, 

(a+ ' a_)W ' \ y (a+i0 -j(t)/Var 1 /a(y (a+iB _ ) (l)), a+ ^ a", 

where A and F( + l0 -) are defined at p. 31) and (14. ip . respectively. Accordingly, the process A( a + a -\ is defined by the 
following harmonizable representation: 



*(«,+,.-)(*) = *(#(*)) / | , gm+1 U [a+ . a - ) (H(t),x)W(dx), te(0,l), (4.3) 
where K(H) is defined in (|1.3j) and 

{1, a+=a~^0, 
Q+e -isign (3; )(g+i) f + Q - c isign( a )(g+i) f ^ + ^_ (4 4) 

((a+) 2 + (a-) 2 - 2a+a- sin(7riJ)) 1/2 

Note that the function H M( a +,a-)(H, x ) is analytic on (0, 1) for all cc € R and (a + , a - ) € R 2 \{(0, 0)}. This definition 
of -X"( + i0 -) induces that for all * € (0, 1) and (a+, a") € R 2 \ {(0, 0)} 

EV,«-)(*) = and EX ( 2 a+,«-)(*) = * 2 " W - 
The following Assumption (C),, is crucial for regularity properties of the MBM, see pQ, [3J. 
(C)„ There exists ry > such that H(-) E C r '{0, 1) and 

0< inf H(t)< sup < min(l,r?). 

te(o,i) te(o,i) 

It is known that Assumption (C) ?) guarantees that the MBM X — A( 11 ) of (| 1 . 3|) is locally asymptotically self-similar 
at each point t € (0, 1) having a FBM Bjjit) as its tangent process at t [T3J and its pointwise Holder exponent coincides 
with H(t) [3J. A particular case (a = (1, 0, • ■ ■ , 0, —1) and m = 1) of the following proposition shows that a FBM B H ^ 
is also a tangent process at t for Xi a +, a -) defined in (14.31) . see Remark [1] and Property [TJ 
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Proposition 1 Let H(-) satisfy Assumption (C) r) . Then, for any a G A(m, q), with m > 1, (a+, or) £ 1 2 \ {(0, 0)} 
and < e < 1/2, £/iere exists a constant Ci(e) > smc/i that for n > 1, 



max 

[ne]<fc<[(l-e)n] 



Var^B^ijC-)) 



The next proposition provides an expansion of the covariance of the process V£ X( a +,a-)- For ease of writing, we 
consider here the case a — a* , m — 2 only. The case of general filter a G A(m, q) is discussed in Lemma |4] (see Section 

ED. 

Proposition 2 Let H(-) satisfy Assumption (C)^ and (a + ,a~) G K 2 \ {(0,0)}. T/ien /or any < e < 1/2, n G N* , 
[ns] <k<k'<[(l- e)n], k' - k > 2q, 



Cov(K I *X (a+ia - ) (-),K l *^(a + , a -)(-)) = A(-,-)vfH(-)vfH(-) 
\ v ' n n / n n n n 



k k' \k - k'\ H (n)+ H ( ! k)- i 
+ B 2 (-,- Y K *' + X n (k,k% (4.6) 

w/iere (recaZZ,) V a * ) = H{^-) - 2H(^)+H(^) are the second-order variations of H(-) and A(t,t'), B 2 (t,t') are 
defined in \1.11^TT^ . The remainder term X n (k,k') in J^.6] ) satisfies the following bound: for any S > there exist 
no G N*, fco G N* such that for any n > n , \k — k'\ > k Q , [ne] < k,k' < [(1 — e)n] 

/ 1 \k' -k\ H ^+ H ^~ 4 \ 

< *(^+ ' n 4 )+H(¥) )■ (4-7) 

Remark 2 From definitions (|7.1HI7.12"]) it immediately follows that A(t, f') and B 2 {t,t') are bounded on (0, 1) x (0, 1) 
and have finite and generally non- vanishing limits 

1 3 

lim A(t,t') = t 2H ^ log 2 i and lim B 2 (t,t') = -- ]J(2H{t) - £). 

Thus, Proposition^ and Lemma|4]allow to obtain lower bounds in (|4.6|) and a decorrelation rate of generalized variations 
of the MBM, sec Corollary [3j (|4.8[) . which is close to optimal. In turn, the decorrelation rate of generalized variations 
determines the decay rate of the variance and the 4th moment of the estimators written in terms of sums of nonlinear 
functions of Hermite rank 2 of normalized generalized variations (see the proof of Theorem [2] also (|6 . 10[) . (1 6 . 11 [) ) . 
Proposition^ and Lemma 0] also permit to construct explicit counter-examples to some earlier results pertaining to this 
issue. In particular, we can show that the expansions obtained in ([15], Lemma 1) and ([11]. Lemma 2) are erroneous. 
A detailed discussion of the above counter-examples is given in the extended version of this paper [B] . Let us note also 
that the Erratum in [16] does not concern ([15], Lemma 1) but another error in the last paper. 

Remark 3 Let us note that in 19 and |7], some correct bounds of the decorrelation rate of VJ? X(-) are obtained, 
but their bounds are less sharp than the bound in (14.8[) . Another correct bound was also obtained in the preprint ([2], 
Lemma 2.4). 

Propositions [T] and [2] help to verify Assumptions (A) K and (B) a for the MBM in (14.31) . 

Corollary 3 Let the conditions of Proposition^ hold. Then for any < e < 1/2, there exist C(e) > and no(e) G N* 
such that for any n > no(s) and any [ne] < k < k' < [(1 — s)n], 

Cor(VfX (a+a - ) (-),VfX (a+a -,(-)) < C(e)( l — — + — L p-Y (4.8) 

V n (a ,a n (a ,a ) \ ^ )) - \ I ^ _ fc | 4 __ fr( |)_ H (^.) „2(^A2) -H ( £) -H{ ) J V ' 

Moreover, 
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(i) ^( +,a-) satisfies Assumption (A) K with to = 2 and any < k < min {(77 A 1) , 2(77 — sup tg ( ^ 

(m/ ^( a +,a-) satisfies Assumption (B) Q wrf/i to = 2 and any < 8 < 7/2, 2 > 7 > 1/2 satisfying the following 
condition: 

7-0 < minf 2 ^ 2 )" 2 ^ , 4-2ff(f)), /or any i G (0, 1). (4.9) 

v 1 — a y 



Remark 4 The presence of the second terms on the r.h.s. of (|4.5I) and (j4.8[) indicates that dependence properties of 
increments of the MBM arc quite sensitive to the smoothness parameter 77 of the Hurst function. These terms have 
a negative effect on estimation of H(-) and the convergence rates, by imposing restrictions on the bandwidth a, see 
Proposition [3] below. [23] argued that the dependence properties of the MBM are rather peculiar and proposed a 
different class of multifractional Gaussian processes defined via nonhomogeneous fractional integration of white noise. 
These processes are locally self-similar in the sense of (11.11) under less restrictive conditions on H(-) and have better 
decorrelation properties than the MBM [53] . On the other hand, the covariance function of the above-mentioned time- 
varying fractionally integrated processes is not a local function of H (t) and H{t') as in the case of the MBM (see (|7.1|l ) 
and its study is more difficult. Extending Propositions [1] and [2] to time- varying fractionally integrated processes of |23] 
is an interesting open problem. 

Proposition 3 Let X (o+i0 -) be the MBM of (l4~3)) with (a + ,a~ ) £ K 2 \ {(0,0)} and H(-) satisfying Assumption 
(C) r) , 77 > 0. Then, with a = a* and £ = IR, IR2, QV or QV2, 

(i) For any t £ (0, 1) and any a £ (0, 1), H^l(t) -A H(t). 

(ii) For any t £ (0,1), H^l(t) ^> H(t), provided H(t) < rj - § and max(0, 1 - 4((?7 A 2) - H(t))) < a < 1/2 hold. 

n— >oo 

(Hi) For any e > 0, sup £<t<1 _ c \H^i a{t) — H{t)\ — Oj^n^ 2 ^ 30 - 1 / 4 ), provided a and r\ satisfy 

a > sup (l-4((r; A 2) -#(<))), (4.10) 
te(o,i) 

1 2 2 

sup Hit) < 77- — and < a < -. (4.11) 

te(o,i) 12 3 + 4(77 A 2) 3 

(iv) For any e > and S > 0, sup £<t<1 _ £ \Hn £ l(t) — H(t)\ = 0(n~( 1 ~ 3a ~' 5 - ) / 4 ) a.s., provided a and rj satisfy \^.1Q^ , 

sup H{t)<ri-\ and ] < a <\. (4.12) 

{e(o,i) 6 3 + 4(t?A2) 3 

(v) For any u G W and i x , • • ■ ,t u G (0, 1), t 4 7^ ^ (i 7^ j) the multidimensional CLTs p~J7| ), fUfJ), iTO]) and ifgJTJj) 
Zio/d, provided the matrix-valued functions r satisfy the condition of Corollary [JJ and £/ie parameters a and 77 
satisfy and 

a > max{-— 1— -, 1-4((77A2)- sup (*))}. (4.13) 
ll + 2(??A2) 4e(0jl) y J 



Remark 5 Roughly speaking, conditions in (ii) - (v) require that 7/ and H (t) are sufficiently separated, or the difference 
7/ — H (t) is large enough. Else, it may happen that a satisfying these conditions does not exist. However, if 7/ > 
(3 + \/4T)/8 ~ 1.175 then a satisfying (ii) - (iv) exists for any H(-) G C(0, 1) with sup tg ( n i/(<) < 1, while condition 
(|4.13j) in (v) reduces to a > 1+2 ^ A2 ) ■ 

Remark 6 In the case of estimators £ = QV^ or QV2, Proposition [3] can be easily extended to a general filter a £ 
*4(m, q), m > 1. For m = 1, the corresponding results in (ii) - (v) can be proved under weaker conditions on a, 77 but 
the convergence rates are worse. On the other hand, if a £ A(m, q), m > 3, the corresponding conditions on a, 77 in (ii) 
- (iv) are the same as in the case a = a* except that ?yA2 in (|4.10[) and (|4. 1 2[) can be replaced by 77 Am. In (v), condition 
(|4.13p holds. Note also that this proposition also holds for a process (<r(i) ^( a +. a -)(' : ))tG(o,i) when a £ C v (0, 1). 
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Remark 7 For a class of harmonizable Levy processes, Lacaux ([T5], Theorem 4.1) obtained the weak consistency 
rate of the QV estimator H^a\t) — Hit) = O p {n~ p \ogn) with p < r/ — H(t),p < in the case when H(-) £ 

C"'(0, 1), < n < 1. For the MBM the above result follows from the CLT in Proposition [3] (v) with p = ±^ satisfying 
p < 2(rj — H (t)) and p < . This means that concerning the weak consistency of the point- wise QV estimator of 
H(t), our results for the MBM are more accurate than those in [19]. Moreover, the last paper does not discuss strong 
and uniform consistency rates and the CLT as in Proposition G3 

Remark 8 Bcgyn [8] obtained strong consistency and asymptotic normality of sums of quadratic variations for general 
non-stationary Gaussian processes under different assumptions on the covariance function which exclude the case of 
MBM. In particular, if r\ < 2 then the covariance R(t,t') = EX(t)X(t') in (|7.1[) is not twice differentiable in t or t' 
and the derivative d 4 R(t,t')/d 2 tdt' 2 for t ^ t' does not exist, contrary to what is assumed in ([8 , Theorem 1 and 
2). In the case of quadratic variations localized in a neighborhood V n>a (t) of given point t £ (0, 1) as in the estimator 
Hnja if) (]2.7p . the imposed conditions in ([8], eqs. (5), (21), (23)) exclude the appearance of the first term of the order 
0(n _2 ' 1,A2 - ) ) on the r.h.s. of (|4.6j) which cannot be ignored and which plays a crucial role in our asymptotic results. 

Remark 9 Since a depends on -q and H (t) that are not available, it is a problem to apply the CLTs of Proposition 
[3] in concrete situations. However, when rj > 2 it is possible to select any a > 1/5 and even to use n 1 / 5 log 2 n instead 
of n a ; then the convergence rate of the CLT is n 2 / 5 up to a logarithm term. Note that [TS] proposed a procedure by 
minimization of the MISE to select an optimal a and then an adaptive estimator of H(t) is H^^\t). 



5 Simulations 



All the softwares used in this Section are available with a free access on http : //samm . univ-parisl . f r/- Jean-Marc-Bardet 
(in Matlab language). 

We use the original version of the MBM X in f| 1 . 3[) . Since trajectories of the MBM in our simulation study are generated 
(for a given Hurst function H(-)) using the Choleski decomposition of the covariance matrix, the number of observation 
points limited to n = 6000. Although this data length may appear rather small in the present context, some interesting 
features can nevertheless be noted. Three cases are considered: Case 1: H(-) is a smooth function; Case 2: H(-) is a 
trajectory of an integrated FBM with Hurst parameter < h < 1 and independent of X, therefore H(-) £ C v ~ with 
r\ = 1 + h £ (1, 2); and Case 3: H(-) is a trajectory of a FBM with Hurst parameter < r] < 1, independent of X. 
For each local Hurst function H(-), Monte-Carlo experiments are realized from 100 independent replications of observed 
paths (X(l/n),X(2/n), ■ ■ ■ ,X(1)) of the MBM (TQ|) for the following choices of parameters: 

• a = 0.2, 0.3, 0.4 and 0.5; 

• p = 5 (= the number of dilatations of H^a^ and , — the regression length providing Hn^) in all cases; 

• a = a* = (1, —2, 1) and therefore m = 2 in all cases. 

Each estimator of H(-) is computed for t = {n~ a , + 0.01, • • • , min(l — n~~ a ,n~ a + 0.99)} and therefore an approx- 
imation of VMiSE = ( / E(H n (t) - H(t)) 2 dt) 1 can be computed. Here there are the results of simulations: 

Case 1: H(-) £ C°° 

We have chosen 4 different functions 1 1— » H(t). These functions are: 

• Hi(t) = 0.6 for any t £ (0,1). 

• H 2 (t) = 0.1 + 0.8t for any t £ (0, 1). 

• H 3 (t) = 0.5 + 0.4sin(5r) for any t £ (0, 1). 



Hi(t) = 0.1 + 0.8(1 - t) suT(lOt) for any t £ (0, 1 
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From -ffi(-) to H^{-), the functions H(-) display more and more ample fluctuations and less regularity, even if these 
function are all C°°(0, 1) functions. Figure [T] provides two kinds of graphs (the mean trajectory and single trajectory) 
of the two estimators H^a^ 1 and of the function H4 for n = 6000 and two values of a (a = 0.3 and 0.4). The 

corresponding MISE are given in Table [5] 



Case 2: H(-) G C~, 1 < 77 < 2 

We have chosen 77 = 1.5 but other simulations with different values of rj lead to similar conclusions. To ensure impartial 
results, we chose to compute the estimators for 50 different functions H(-) generated from 50 trajectories of integrated 
FBM. Then 100 independent replications of processes generated with each function H(-). Table [3] contains the empirical 
MISE computed from these 50 x 100 processes. Figure [5] provides the trajectories of Hnfa (t) an d Hn,a 2 \t) for one 
of the 50 differentiable Hurst functions H (•) G C l l5 ~ for n = 6000 and a = 0.3, 0.4. 



Case 3: H(-) e C ? '~, < 77 < 1 

We have chosen here r\ = 0.6. As in the case 77 = 1.5~, we generated 50 different trajectories H(-) obtained from 
trajectories of FBM with parameter H = 0.6. Then for each function H(-), 100 independent replications of MBM with 
Hurst function H(-) are generated. Table [5] contains the empirical MISE computed from these 50 x 100 processes. An 
example of a graph of a function H(-) e C 6 ~ and the mean trajectories (obtained from the 100 replications) of both 
the estimators are drawn in Figure [3] 



Conclusions of simulations: 

1. In any cases of function H(-), for n = 2000 or n = 6000 and a = 0.2, 0.3, 0.4 or 0.5, the empirical MISE of 
Hn^a 2 \t) and Hn^ 2 \t) are smaller than the ones of Hi^a (t) and Hnfx(t), respectively. This theoretically 
corresponds to the Gauss-Markov Theorem which is also empirically satisfied. Note that PGLSE (t) allows 
for certain cases to divide by 2 the empirical V MISE of Hn$ (t) , while the gain obtained with H^^f 2 ^ (t) instead 
of Hifja(t) is more limited (less than 10%). 

2. In agreement with the theory (see the remark at the end of Section^) in the case H(-) € C v with 77 = 1.5~ or 77 = 00, 
our simulation suggests to select a greater than (1 + 2(77 A 2)) _1 for "optimal" results. More precisely, it seems 
that a good choice is a ~ 0.3 for Hn^} and Hn^\ and a ~ 0.4 for and H^a 2 ^- This empirical rule seems 
also valid for H(-) € C 6 even if the theoretical choice should be a > max{5/ll , 2(2sup t6 ( ^ H{t) — 0.7)}. 
However a bias clearly appears for a = 0.3 when H(t) is close to 77, while this bias is reduced when a = 0.4 
(see Figure [3]). The tables also confirm that, with the above choice of the bandwidth a, the accuracy of both 
estimators increases according to the theoretical convergence rate n^ 1- "^ 2 (~ 71 1 / 3 ) as n increases from 77 = 2000 
to 77 = 6000 (for instance for a = 0.3 the empirical V MISE of Hn^ 2 \-) seems to be divided by 1.5 ~ 3 ' 35 ). 



3. Mutatis mutandis, the values of the empirical \J MISE are quiet the same for any functions H(-) £ C with 77 > 1. 
However, even if this observation seems to be able to be extended for 77 = 0.6, this does not hold since in this case 
the function if (•) is required to vary from to 0.6 while for 77 > 1 it varies from to 1. 

4. Finally, in all cases, the best estimator is clearly Hn,a with the choice a = 0.3 (even when if is a constant, the 
empirical V MISE of for a = 0.3 is smaller than the one of H^a 2 ^ for a — 0.4). However the asymptotic 
variance of Hn^ 2 \t) is a little larger than the one of H^f 2 \ But, as it can be seen in Figure [T] and Figure EJ 
Un,a 2 ^ is nearly unbiased when a = 0.3 (and sometimes when a — 0.2) while the one of H^f 2 " 1 is still biased for 
a = 0.3. Then Hna 2 ^ can be used when a = 0.3 with an asymptotic variance varying with n^ 1 "™)/ 2 = n~ - 35 
while Hk Q r ] requires a = 0.4 and an asymptotic variance varying with 71 • . 
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Hi(t) = 0.6 


a 


0.2 


0.3 


0.4 


0.5 


n = 2000 


V MISE for Hnfa 


0.044 


0.055 


0.073 


0.104 




VM/SB for H£? a ' 


0.041 


0.051 


0.069 


0.096 




V MISE for JiA.a 


0.111 


0.137 


0.186 


0.260 




. / 71 /TTCITTl £ 7-7- (Z/?2 ) 

V MISE for ii„, Q ' 


0.061 


0.077 


0.106 


0.145 


n = 6000 


VM/SB for • 


0.026 


0.035 


0.053 


0.079 




V M/SB for HL Q / 2) 


0.025 


0.033 


0.050 


0.074 




VmISE for M!q } 


0.065 


0.091 


0.128 


0.202 




V M/SB for Hum 


0.037 


0.049 


0.076 


0.115 


H 2 (t) = 0.1 + 0.8* 


a 


0.2 


0.3 


0.4 


0.5 


n = 2000 


/- ^-/nT/ 1 ! 

V MISE for i^ Q Q y ) 


0.170 


0.076 


0.075 


0.101 




VM/SB for HZg 2) 


0.170 


0.073 


0.072 


0.096 




V MZSB for JiA.a 


0.115 


0.143 


0.184 


0.247 




V MISE tor ii A, a 


0.059 


0.071 


0.098 


0.135 


n = 6000 


VM/SB for r ) 


0.115 


0.045 


0.051 


0.074 




VM/SB for iJ^a ' 


0.114 


0.044 


0.048 


0.070 




VMISE for C' 


0.070 


0.094 


0.134 


0.195 




V MISE for iiA.a 


0.036 


0.046 


0.069 


0.103 


H 3 (t) = 0.5 + 0.4sin(5i) 




0.2 


0.3 


0.4 


0.5 


n = 2000 


VM/SB for > 


0.362 


0.125 


0.084 


0.102 




VM/SB for ! 


0.362 


0.123 


0.080 


0.096 




V MZSB for JiA.a 


0.129 


0.133 


0.171 


0.229 




V M/SB for Hn.a 


0.093 


0.071 


0.091 


0.124 


n = 6000 


VM/SB for M Q a ^ 


0.260 


0.078 


0.056 


0.077 




VMISE for Jf4,Q - 


0.260 


0.077 


0.052 


0.072 




VMISE for M!a } 


0.078 


0.089 


0.125 


0.180 




V M/SB for ^a 2) 


0.057 


0.047 


0.065 


0.097 


H 4 (t) = 0.1 + 0.8(1 - t) sin 2 (10t) 


a 


0.2 


0.3 


0.4 


0.5 


n = 2000 


/- — — — - 
VM/SB for ' 


0.321 


0.165 


0.121 


0.120 




VMISE for i/4?a ' 


0.320 


0.164 


0.117 


0.112 




V MISE for Mfa" 1 


0.178 


0.138 


0.160 


0.210 




VMZSB for H^a 2) 


0.165 


0.098 


0.091 


0.112 


n = 6000 


VMISE for ^ Q Q y) 


0.251 


0.136 


0.074 


0.084 




Vm/SB for HL Q / 2) 


0.251 


0.135 


0.071 


0.078 




VmISE for fl*™* 


0.158 


0.088 


0.115 


0.164 




VMISE for £<S 2) 


0.148 


0.062 


0.067 


0.091 



Table 1: Values of the (empirical) MISE for estimators H^a \ H^a 2 \ Hnf* and Hn.a 2 ^ of C°°-Hurst functions for 
n e {2000, 6000}, a e {0.2, 0.3, 0.4, 0.5} and a = (1, -2, 1). 
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alpha=0.3 



alpha=0.4 




0.4 



0.2 0.4 0.6 Q.t 

t 




alpha=0.3 



alpha=0.4 




0.2 0.4 0.6 O.f 

t 




Figure 1: Estimates of the function Hiit) = 0.1 + 0.8(1 — t) sin 2 (10i)) with t € (0, 1) for n = 6000 and a = 0.3 and 
0.4 (from left to right). The top row represents the mean trajectories of Hnj£f 2 \t) and Hn,a (t) obtained from 100 
independent replications of MBM with the above function H(-). The bottom row represents a trajectory of H^^ 2 \t) 

' — ' (I R2) ~~~~ (OV2^i 

and Hn.a (t) obtained from a trajectory of a MBM with the same function H(-). The graphs of H(t), H)Ca (t), o,nd 



H, 



(IR2) 



(t) are in black, blue and red, respectively. 
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A 1 OA 

0.139 


A 1 /I 1 
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A 1 7CC 

0.175 


A OO O 
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V MISE for SL 1 ^ 
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0.077 


0.093 


0.128 
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V MT^F for H^ IR2 ^ 


u.UOO 


U.U ( o 


u.UJU 


U. ±o<j 


n = 6000 


for ffj$» V) 


0.129 


0.067 


0.057 


0.074 




VMKE for M Q Q y2) 


0.130 


0.067 


0.056 


0.071 




for 


0.096 


0.106 


0.143 


0.201 




for M™ 2) 


0.066 


0.052 


0.067 


0.103 



Table 2: Values of the (empirical) MISE for estimators H^a \ Hn^a 2 \ Hn^} and of C -Hurst functions for 

r] = 1.5 (up) and rj = 0.6 (down), n € {2000, 6000}, a £ {0.2, 0.3, 0.4, 0.5} and a = (1, -2, 1). 




Figure 2: Trajectories of Hn,a (t) and Hn, a (t) for one of the 50 differ entiable Hurst functions H(-) G C 1,5 ~ /or 
n = 6000 anc? a = 0.3, 0.4 (from left to right). The graphs of H(t), Hr^ 2 \t) and Hn^\t) are respectively in black, 
blue and red. 
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Figure 3: Mean trajectory (from 100 independent replications) of H^f\t) and Hn$(t) for one of the 50 differ entiahle 
Hurst functions H() G C°' 6 ~ for n = 6000 and a = 0.3, 0.4 (from left to right). The graphs of H(t), H^Sof 2 \t) and 
ifnfa (t) are respectively in black, blue and red. 

6 Proofs of Theorems Q], [2] and Corollaries [2] 

The proof of Theorem Q] uses the moment inequality in Lemma [IJ below, which is a particular case of a more general 
moment bound in [4]. 

Let (Y"i,--- , Yjv) be a collection of Gaussian vectors Yt = (1^ , • • • , Y/" ) G with zero mean EY f = and 
non-degenerated covariance matrices S t = (Cov(Y^ u \ Y^)) < , having a joint Gaussian distribution in K"^. Let 

e G [0, 1] be a fixed number. Call (Yl, • • • , Yw) s— correlated if | Cor (l^ , " ') | < £ for any t ^ s,l < t, s < N and 
any 1 < u, v < v. Finally, denotes the sum over all distinct integers 1 < ti, ■ ■ ■ ,t p < N, U ^ tj (i ^ j). 

Let \x\ = (E J =i( a;( ' ) ) 2 ) 1/2 denote the Euclidean norm in M", A = (dy) a real v x is matrix, A T the transposed 
matrix, / the unit matrix, and ||A|| := supi^i—j \Ax\ the matrix spectral norm, respectively. Lemma Q] is a particular 
case of the moment bound established in Bardet and Surgailis (2011, Corollary 1). 

Lemma 1 Let (Yi, • • • , Yjv) € M." N be an e— correlated Gaussian vector such that 

max IIS^ 1 !! < c max (6.1) 

l<t<JV 

for some constant c max > 0. Let Gjj : — > R, 1 < j < p (p > 2), 1 < t < N be measurable functions such that 
||Gj,t|| := E|Gj i t(l r f)| 2 < co. Assume that for some integer < a < p, 



and that £z^ 2 c max < ~ i ■ Then 



E[G jtt (Y t )] = 0, E[Y t G jtt (Y t )] =0, for all 1 < j < a, 1 < t < N (6.2) 
i 

up—I ' 

\nG 1 , tl (Y tl )---G p , tp (Y tp )}\ < C{e,p,a,c max )KN^Q%, (6.3) 

where 

Qn := max £ i max \EY^Y 8 ^\ 2 , (6.4) 

l<t<N z — ' l<u,v<v 
l<s<N,s^t - ' - 

K := n?=i maxi< t <jv ||Gj,t|| , and the constant C(e,p, a, c max ) depends on e,p, a, c max only. 
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In what follows, we use Lemma [T] with 2 < p < 4, f = 2 only. Let 

Q (A 1 X- m \VfZ(k/n) + VfZ((k + l)/n)\ 

»»,«(*) ke £t (t ) l y « z ( fc / n )l + IK? ^(0 + i)/™)! 

Thus, ff^ ) W=A^ 1 (5 n (t)). 

Lemma 2 Let Z = (Z(t))te(o,i) &e a zero-mean Gaussian process satisfying (A) K wi£/i m = 2. TTien /or any < e < 
1/2, 

sup |E£ n (i) - A 3 (H(t))\ = o(n- K ) + 0(n- 1 +n- Q (" A2 )). (6.6) 

e<t<l-e 

Proof of Lemma fj[ The first relation in (|6.6[) follows from 

max |E»fc(fc)-A 2 (£r(fc/n))| = ( n -»), (6.7) 

lne/2J<fc<[(l — £/2)n] 

sup Wn,a(<)| E {A 2 (ff(fc/n))-A 2 (if(i))}| < Cn"^ A2 ) (6.8) 
feev„, a (t) 



e<t<l-e 



To show flCTT]), write Ery n (fc) = E[?A(K l Q *Z(fc/n) , Z((k + l)/n))] = A(Cor (F n °* Z(fc/n) , V n a *Z((fc + l)/n)) where 
A(p) := ~ arccos(-p) + ^y^^^log (j^), see ([5], p.777), also (|2~5|) . From Property!!] 

max n"|Cor(K°*Z(A/n), V; *Z((fc + l)/n)) - pa(H(*;/n))| — >0. 

[ne/2]<fc<[(l-e/2)n] 1 ' 1 n->oo 

Since K.2(H) = K(p2{H)) and the function p G (—1, 1) n- A(p) is analytic we deduce (|6.7I) . 
To show (|6.8I) . write 



A 2 (H(k/n))-A 2 (H(t)) = A' 2 (H(t))(H(k/n)-H(t))+0((H(k/n)-H(t)) 2 ) 

= A' 2 (H(t))(H(k/n) H(t)) + G>((£ - *) 2( " A1) ) 
= A' 2 (H(t))(H(k/n) - H(t)) + 0(n~ 2a( ^) 



for fc G V n , a (t). Hence 

E {A 2 (ff(A/n))-A 2 (ff(*))}| < C«-i(t)| E {H(k/n)-H(t)}\+Cn- 2 ^D 
fc6V„,a(t) feev„, a (t) 

< C ( n -a(,Aa) +n -l +n -2a(l J Al)) > 

where for 77 > 1 we used H(k/n)-H(t) = H'{t){{k/n)-t) + 0{n- a ^ A ^) and £ fceVn a (t) {(&/»*)-*} = O^ 1 ). 

This proves f|6.8|) and Lemma [21 too. □ 

Proof of Theorem^]] (i) Recall Hn,a (t) = A 2 _1 (5 I1 (t)). Note that the map H i->- A2{H) is continuous and strictly 
increasing on [0, 1] and the inverse map s H> Aj (s) is continuously difierentiable on s G [A2(0), A2(l)]. Therefore, since 
and S n [t) G [0, 1] a.s., so part (i) follows from S n [t) — > A2(H(t)) which is a consequence of 

n— too 

ES n (t) — ■» A 2 (JT(t)) and S„(t)-ES„(t) A 0.. (6.9) 

The first term of (|6.9p is obtained from Lemma [21 The second member of Lemma [5] can be deduced from the following 
bound: 

Q 

E(S n (t)-ES n (t)) 4 < Cn- p , for some p := 2(1 - a)(l - -) > 0, (6.10) 

7 
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from Assumption (B) Q with a e (0, 1). 

We shall prove a slightly stronger bound: there exists a constant C < oo such that for any integers < s < t < 
n, \t-s\< 2n 1 - a 



E( Vnik)) 4 < Cfi-sfn^, where fj n {k) := r] n (k) - Er) n (k). (6.11) 

k=s+l 



Relation (|6~TTj) implies (joUOt with p := 4(1 - a) - 2(1 - a) - 2e( ! y ' a) = 2(1 - a)(l - |). 

In order to apply Lemma [T] to prove ()6 . 1 1|) . we rewrite the ratios rj n {k) = ip{Y n {k)) as functions of Gaussian vectors 
Y n (k) = (r„ (1) (fc),yj 2) (fc)) e K 2 , where 



(1) V: Z(k/n) (2) _ V:Z((k + l)/n) 

K " {k) - a n {k) ' *» - <7„(fe + l) • [t ' 

a 2 n (k):=V a r(V 7 fZ(k/n)). (6.13) 



Also define 



U V a *B H(k/n) (0) V a *B H{k/n) (l) 

y/Vax(V°*B H(k/n) (0)y ^{V a *B H[k/n) {l)) 

Note V( - 2/ (1) ,~l/ 2) ) = ^(y (1) ,y (2) ), implying E[Y„(fc)V'(K(fc))] = 0; c.f. condition of Lemma [TJ Property 

[T] implies that the covariance matrices T, n (k) of the Y n (k)'s approach as n —> oo uniformly in [ne] < k < [(1 — e)n] 

the corresponding non-degenerate covariance matrices S„(fc) := ( P2( ( / )) ] f Gaussian vectors 

' \ P2(H(k/n)) 1 J 

W n (k) = (Wn\k), W n 2 \k)^ and therefore these matrices satisfy the bound in (|6.ip . viz., 

max n K ||E„(fc)-S„(fc)|| -> 0, max ||S„(fc)- 1 || < c max (6.15) 

l<fc<ra-2 l<fc<n-2 

with some c ma x < oo independent of n. 

Next, we "decimate" the sum on the l.h.s. of (|6. 1 1|) so that the "remaining" vectors Y n (k)'s are e'— correlated, as 
follows. Let I = l n be the sequence of integers increasing to oo at a rate o(n 1_Q ) which is specified in (|6.21|) below. 
Write 

t £-1 

R n (s,t) = $3^0"). ^(i) : = 2 ^" (fc) - (6 ' 16) 

fc=s+l j'=0 s<k<t: k=j(mod t) 



Then 



where 



E(i?„(M)) 4 <** max E(i? M 0')) 4 and E(i?„.,(j)) 4 < C fe 4 + E 3 + £2 i + £2 2) . (6-17) 



£4 : = X! | E [*7n(fcl) ' ' ■^(^4)] |, £ 3 := X! | E ['7n(*l)'/n(* ; 2)^(A;3)]|i 

£2,1 := £'l E [^fo)*M*3)]|, £2,2 := EE^ft^fe)], 

where (respectively, J^) stands for the sum over all different integers s < k p < t, k p — u (mod £), k p ^ k q (p q) 
(respectively, over all s < k\, ki < t, k\ = u (mod £), k^ — u (mod £)). But, since |?7„(A:)| < 1 for all k, n, 

E2,i + E 2 ,2<C(\t- S \/£) 2 . (6.18) 

Let us estimate ^ 3 and J^ 4 . Let 

QnAl) ■= , „ max v J4 , S Pn{hi,k 2 ), (6.19) 

s<&2 <t,k2— j(mod £),k2^ki 
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where 

p n (h,k 2 ):= max \E[Y^ (k^Y^ {k 2 )\ \ 

and r„ (1) (fc),r„ (2) (A;) are denned in (|6.12[) . Assumptions (B) Q and (A)o imply the following bound. There exists a 
constant C\ < oo such that for any n > 1 and any [ne] < k\ < k 2 < [(1 — e)n] 

Pn(ki,k 2 ) < ^n^-^dh-k^An 1 - )^- (6.20) 

see also [4]. Now, given C\ > 0, define 

C 2 :=d/e' and £ := [Qnt 9 ^" 1 -"']. (6.21) 

Relations (|6.20[) and (|6.21l) imply that for any < j < £, the Gaussian vector (Y n (k), k — j (mod £), k = 1, • • • , n) is 
e'— correlated. We choose e' < l/(28c max ), thus guaranteeing condition e' < 1 / '(c max ^ 2 (yp — 1) of Lemma[T]for v = 2 
and p = 3, 4. 

Using ([6720]) . (|OT|) and j>l/2,\t- s\ < 2n 1 ~ a we obtain 

max O n t(i) < C ( — ) max > \ki — k 2 \ 7 



9(l-a) 



OC 



< 2C(^^) 2 ]>>'r 27 < C (6.22) 
fc'=i 

Thus, the sums £ 4 and £ 3 can be bounded using Lemma Q] with N = [\t — s\/£]. Accordingly, for £ 4 with p = a = 4 
we have 

E 4 < C (|* - ( max Q M (j)) 1 < C (|t - s|/£) 2 . (6.23) 

Similarly, for with p — 3 and = 2 we get 

£ 3 < C (\t - s\/£f-i ( max Q M (j)) 1 < C (\t - S |/£) 2 . (6.24) 

Relation (|6TT|) follows by combining (foTTTj) . (|6l8|) . (j03]) and (f6T24|) . 

(ii) Similarly as in the proof (i), it suffices to show 

S n (t)-ES n (t) ^ 0. (6.25) 

n— >oo 

In turn, (|6.25[) is a consequence of (|6 . 10[) with p > 1 (since a < 2 J-y-8) ) and Chebyshev's inequality imply that for any 
< S < (p- l)/4 

OO OO OO 

^P(|5„(t)-E5„(i)| >rc" 5 ) < £ 7i 45 E(5„(t) - ES n (i)) 4 < C^n 4 ^" < oo. 

n— 1 n— 1 n— 1 

Hence, (|6.25j) follows by the Borel-Cantelli lemma. 

(iii) We first estimate sup e<t<1 _ e \S n (t)) - A 2 (H(t))\, where S n (t) is defined in Split S n (t) - A 2 (H(t)) = {S n {t) - 
ES n (t)) + (ES n (t) — A 2 (H(t))). Since u„ = w„, Q (t) ~ 2n 1_Q do not depend onte (e, 1 — e) for n large enough, so 

sup \S n (t) - ES n (t)\ < v~ x max |i?„(s, s + v n )\ < v" 1 V] max |i?„(s, s + u„)|, (6.26) 

e<i<l — e 0<s<n— ^™ z — * (i—l)v n <s<iv n 

l<i<n/v n 

with R n (s,t) as defined in (|6.16p . Therefore for any a > 

P( sup IS^i) - E£ n (i)| > a) < V P( max \R n {s, s + v n )\ > av n ) . (6.27) 
e<t<l-£ , ^r~^ , (i-l)v n <s<iv n 
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Since \R n (s,s + v n )\ = \R n ((i - l)v n ,s + v n ) - R n ((i - l)v n ,s)\ < \R n ((i - l)v n ,s + v n )\ + \R n ((i - l)v n ,s)\, therefore 
P( max \Rn(s, s + v n )\ > av n ) < 2P( max \R n ((i - l)v n , s)\ > av n /2). 

{l—\)v n <~s<iv n (i—l)v n <~s<(i-\-\)vn 

We will use the following fact from Billingsley ( [14] . Theorem 12.1, see also (12.5), (12.10)). Let £j, 1 < j < m be 
arbitrary r.v. 's. Suppose there exist nonnegative numbers Uj > 0, 1 < j < m such that 

s s 

E( ^ ( E "if' 0<i<.s<m. (6.28) 

j=t+l j=t+l 



Then for any e > 

s mm 

p ( ^ m I E ^ e ) ^ K ^ ( E u i) 2 + p ( I E & \ ^ e / 2 ) > ( 6 - 29 ) 



where K > is an absolute constant. 
We shall apply (|6T29l) with 



1 /2 At 1 -") 

£j := Vn((i ~ +j)> m 3 := C 7 n -» , m := u„ - 2n 



l-a 



where C > is the same as in (I6.11[) . The validity of condition (|6.28|) follows from the bound (|6.11[) . Hence from (|6.29[) . 
(|6.11[) and Chebyshev's inequality we obtain 

P( max \R n {{i-l)v n ,s)\> av n /2) < C(aw„/4)~ 4 ^n M V^ i + P(|i?„((i - l)v n , iv n )\ > av n /A) 

(i—l)v n <s<(i-\-l)v n 

< CaT^rr^ + C(av n )- A V\R n {{i - l)v n ,iv n )\ A 

_„ 26(1- a) 

< Ca v n n 1 

Hence and from (|6.27|) we obtain 

P( sup \S n (t)-ES n (t)\>a) < Ca- 4 n a - 2{1 - a){1 -^\ (6.30) 

e<t<l-c 

where C > does not depend on a and n. Taking a = <5n~^ with jj, given in (|2.16[) results in P(sup e<t<1 _ e |S n (t) ~ 
ES„(t)\ > Sn-^ 1 ) < CVT 4 , V<5 > 0, or 



sup \S n (t)-ES n (t)\ = O p (n~»). (6.31) 

e<t<l-e 



Next, from Lemma © we have 



sup \ES n (t) - A 2 (H(t))\ = o(n- K )+0(n- a ^ 2 )) = 0(n _ "), (6.32) 

e<t<l-e 

since 2a < 1. Finally, (|6.32[) and (|6.31l) imply 

sup \S n {t) - A 2 (H(t))\ = O p (n-"). (6.33) 

e<t<l-e 

The result in (iii) follows from (|6.32|) and the properties of the inverse map s h-> A^" 1 (s) mentioned in the proof of (i). 

(iv) Follows from the Borel-Cantelli lemma, (|6.30[) with a = n^ 1 -^ 1 ^^ and (|6.6p with /i replaced by fix (which in turn 
follows from (|6~7|) . and ([5T5])). Here, we use the facts that 2(1 - a)(l - {6/j)) -a- A(fj, 1 - 8) = 1 + 45 > 1 and 
a(?7 A 2) > fix- Theorem Q] is proved. □ 



The proof of Theorem [2] below, uses the following lemma which is a particular case of a more general result in [4] 
(see also [5], Theorem A.l). 
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Lemma 3 Let Yk,N = (¥k%i " ' i ^few) ^ ^"i 1<&<-^V, -'VGN&ea triangular array of jointly Gaussian vectors, with 
zero mean EY^jv = and non- generated covariance matrices E^at — (Cov(Y fe ' t 2-, ij.^)) 1<u v<v , I < k < N, N e N. 
Assume that there exists a function p : Z — > R smc/i £/ia£ 



V(j,fc)G{l,---,7V} 2 , max EY$Y$ < - fc)| with ]T |p(j)| 2 < oo. (6.34) 



Moreover, assume that for any r 6 (0, 1) any JgN*, 



( y j+[WT],w)_ J < J < J 1 11 ! ./ , ./ • ( 6 - 35 ) 

— J — Tl->00 

where Wj = {W^\ • • • , wf- u) ), jeZ is a stationary Gaussian process taking values in W and such that the covariance 
matrix £ := (Cov(VF (tl) , W^ v) )) 1<u . is non-degenerated: det £ ^ 0. in addition, assume that 

max ||£fciv-£|| -> (N ->• oo). (6.36) 

l<fe<AT ' 

Lei G : — > K 6e a bounded measurable function such that 

E[Y kiN G(Y kiN )] = 0, Vl<fc<iV, TVeN. (6.37) 
JTien, twt/i a 2 := ^ Cov(G(tY ), < oo, 

AT 

Ar- 1 / 2 V){G(n,Jv)-EG(n,Jv)) A .A/(0,a 2 ). (6.38) 

k=l 

Proof of Theorem [H For notational simplicity, we restrict the proof of the p— dimensional CLT in (|2.17[) to the case 
p = 2. Let < ti < t 2 < 1 be arbitrary time moments and S n (t) be defined as in (|6.5p . First, we prove the following 
bivariate CLT for (SV^i), S n (t 2 )): 

v. iffn ^ ( o- 2 (H(h)) x 

^(H(t 2 )) 



V2n^(S n (t 1 )-A 2 (H(t 1 )),S n (t 2 )-A 2 (H(t 2 ))) AjV(Q,r), T := ^ ^ " , (6.39) 



where 271 1 " ~ w n = #V n , a (ti),i = 1,2, for all n large enough. Relation (|6.39[) follows from E5* n (ii) — A 2 (H(ti)) 
^ n (i-a)/2^ ( wn j cn i s a consequence of (I6.6[) and the inequalities for k and a in the statement of Theorem [2]) and 



(2n 1 - a ) 1 / 2 (S n ~ES n ) -A Af{Q,a 2 ), (6.40) 

n— >oo 

S n := c.Sih) + c 2 S n (t 2 ), a 2 := efc 2 ^)) + c\o 2 (H (t 2 )) , (6.41) 



where ci, C2 are arbitrary numbers. 

To show (|6.40p we use Lemma O with N := v n , v = 4. Let ip{y^\ j/ 2 ') = rj^rmq^sjj as above and rewrite 



Sn = N-^GiY^), G(y):=cMy {1) ,y {2) ) + c2*P(y (3 \y {4y ), y=(y^,i 2 \i 3 \y < - 4 y)eR i , 



k=l 

i - ■ ■■= ( Y£% • • • , Y k W) , j ni := K - n l -% i = 1,2, 



yd) , = F (2) . = KfZ(^±l) (3) y°'Z(*±**) (4) . = Vg^(gg^) 



Similarly to (|6. 14|) dchnc 
W k := (wl}\--- ,W^ 4) ) , fceZ, jr* 2 (fl) := Var(l/ a *B ff (0)) = 4 - 4 H , 
(1) ^gwW (2 ) ^VttoM ( 3 ) . = V a * B 2 . mt2) {k) (4) ^X^W 
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where (B it H( t .)(i)),i = 1,2 are independent FBM's with respective parameters Hfa). By definition, (Wk) k£Z is a sta- 
tionary Gaussian process with Wo having a non-degenerated covariance matrix S = (Effi!"' IK ] , E(W ( j"' ) ) 2 = 

V / l<u,v<4 

1, u= I,--- ,4, E[W ( V (2 } ] = ^(^(ii)), E[wfV (4) ] = ^(^(ta)), E[wt } W^] = G, ue {1,2}, v e {3,4}. 
By definition, E(Y fe (i ^) 2 = 1, % = 1,3. From Property HJ it follows that maxi< fe <jv \E(y£%) 2 - 1| 0, i = 2,4 and 



max 

l<k,j<N 



Cov(Y^%Y^)-Cov(W^M a) ) -> (6.42) 



for it, w e {1,2} and u,v € {3,4}. For u e {1,2} and w € {3,4} relation (|6.42|) follows from (B) and (|2.10p . implying 
for the above u, v that maxi<k, j<N 

\Cov(Y$,Y$)\ < C(n\t 2 - ti|)" 7 -> 0. Relation (jQ2|) guarantees the finite- 
dimensional convergence in (|6.35|) for the Gaussian vectors (Yfc,jv) and (Wfe) defined above, including the convergence 
(|6.36[) . The covariance bound in (|6.34l) also follows from (B) since 7 > 1/2. Condition (|6.37p follows from symmetries 
Yk,N l( = — Yk,N and G(y) = G(—y). Finally, the expression for er 2 in (|6.4ip follows by independence of Bi ;H n^ and 
B2,H(t 2 ) and the definitions of G(y) and cr 2 in f)6.38[) . This proves (|6.40p and (|6.39p as well. Then, the clt of Theorem 
[2] (|2.17p for p — 2, together with the expression (|2.18p for the limit variance, follows from (|6.39p by an application of 
the Delta- method. Theorem [5] is proved. □ 

Proof of CorollaryUi By definition, 

m^\t)=j2h{ms{t))m^\t), where h{ H ) ■■= Ij^l^Z ' 

where e, :— (0, • • • , 0, 1, 0, • • • , 0) T <E W is a unit vector. Using the fact that Y2i=i bi(H) = 1, we can decompose 



i-1 



H^ 2 \t)-H = W n (t) + U n (t), (6.43) 

where 

Wn (t) = Y,HH){H^{t)-H), U n (t) := J2HH^\t))-bm)(H^(t)-H) 

i=l i=l 

If view of the assumptions on Y,^(H) : the functions H bi(H), 1 < i < p are continuous on (0,1) implying that 
I C/n (t) I < (maxi^^p |ff™'*(t) — H\)5(\Hr^a (t)) — bi(H)\), where S(x) — o(x), x — > 0. As a consequence, in all cases 
(i) - (iv) the term U n (t) is negligible w.r.t. the term W n (t). From the proof of Theorem 1 it is evident that all statements 
of Theorem 1 hold with replaced by any estimator , 1 < i < P- This proves the corollary. □ 

Proof of Corollary \^ From the decomposition (|6.43l) , 

s/2n 1 -"(HV™\t)-H) = V2ni-«(W„(i) + o p (l)). (6.44) 

Clearly, it suffices to prove the corollary with Hn^ 2 \t) — H replaced by W n (t). In turn, this follows from the (p x 
m)— dimensional CLT 

V^((H^%)-H( tj )) 1<i<p ) A ((W^(t 3 )) ) , (6.45) 

V — — y / l<J<m n— >oo \ — — y / l<j<m 

where (W^^' l (Aj')) 1<i<p ) j = li*" i m are independent Gaussian vectors with zero mean and respective covariance 
matrices Y^(H(tj)).~fhe proof of (jOS]) mimics that of Theorem [2] and we omit the details. □ 

7 Proofs of Section |4] 

This section contains the proofs of Propositions [T] and [2] concerning the asymptotic behavior of (co)variances of gen- 
eralized variations of the general MBM -X"( + )0 -) for (a + ,Gf ")eE 2 \ {(0,0)} with local Hurst function H(-) satisfying 
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Assumption (C) . We discuss the general case of niters a := (a , • • • , a q ) e R 9+1 satisfying condition (|2.1I) for some 
m G N*. 

First introduce the following notations. Using Theorem 4.1 of [35], for f,t' G (0, 1) 



^-)(t)V,.-)( 1 ') = Q(H(t),H{t'),t,t'), 



(7.1) 



where Q is a function of 4 variables € (0, 1), ff, if' € (0, 1) satisfying Q(H,H',t,t') = Q(H',H,t',t) and defined 
for t < H by 



with 



L(H, H')(- sin ( A/3(ff, ff')) (t' log t' - t log t - {t' - t) log(t' - t)) 

+ cos(A/3(ii, if')) * 
i fi u (H, if') + L 22 (H, if') t' H+H ' - L 12 (H, ff') (f - t) H+H ' 



if if + if' = 1 (7.2) 
if if + if' ^ 1, 



Ln(H,H') 
L 22 (H, ff') 



L(H,H') cos (Af3(H,H') - §(ff + if')) 

cos (§ (if + ff')) 
L(if, if') cos (A/3(if, H') + %(H + H')) 

cos(f (if + H')) 



Li 2 {H,H'), 



and where 



K(-) is defined at (|T3|) . Aj3(H,H') 



L(H,H') := 
p{H')-p{H), 



(K{H)K(H')) 



n\V 2 



fv( 



H+H' 



0, a+ = a" ^ 0, 



(7.3) 
(7.4) 

(7.5) 
(7.6) 



where Arg(z) € [0, 27r) is the principal argument of the complex number z 7^ (note that £(if, if) = 1 and 
Lij(H,H',t,t') = L(H, if') for a+ = a" 7^ 0). In particular, for any a ^ 



T'\(+H+H' 1 ±/H+H' 



EX (aia) (t)X (OiO) (0 = EX(i)X(t') = -L(H } H')(t H+H +f 



(t' -t) H+H ') 



t < t'. 



(7.7) 



Note that the functions (if, if') 1— > L{H,H r ), (if, if') i-> Af3(H,H') are analytic on (0, l) 2 and hence (if, if') i-> 
Q(H, H' ,t,t') is also analytic on (0, l) 2 for t' ^ t with possible exception on ff + if ' = 1 where the denominator 
cos (f (if + if')) in (O-dOJl vanishes. To show the analycity of (H,H') i-> Q(H,H',t,t') on if + if' = 1, let 
w := if + if ' - 1 7^ 0, < t < t' < 1 and rewrite Q(H, if', t, t') in (JO) as 



Q(H,H',t,t') = ^L(H,H')cos(Al3(H,H'))(t 1+ " +t' 1+ " -(t' -t) 1+ul ) 



(7.8) 



lL(H,H')sm(A/3(H,H')) 



tan(^) 



— ) - * (— ) + (* - *) (- 



Then, the first term on the r.h.s. of (|7.8[) clearly is analytic in uj at cj = and the same conclusion holds for the second 
term because of the analycity of the functions u 1— > tan ("„./ 2 ) an< ^ w ^ w _1 (i" — 1) (t > 0) in a neighborhood of w = 0. 
It is easy to check that the limit of (I7.8[) as uj —> coincides with the definition of Q(H, ff', t, t') in (|7.2j) for if + if ' = 1. 
As a function of four variables, Q(ff, if', t, t') is analytic on (0, l) 4 except for t = t' but Q is continuous on (0, l) 4 . 
For < t < t' < 1, introduce the 2 x 2-matrix of the second order partial derivatives: 



(R%] H ,{t,t>) R%(t,tf 



(7.9) 
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whose elements are the corresponding partial derivatives of the function Q computed at H := H(t), H' :— H(t'), 

R H,H'(tX)--=^p(H(t),H(t'),t,t'), ••• , R$(t,t'):=^(H(t),H(t>),t,t'). 

In particular, for all < t < t' < 1, 

R (2) (t t , ) = { h HH(t),H(t')) sin {[3(H(t')) - P(H(t))) (f - ty 1 if (t) + H(t') = 1 

*' \ -\L l2 {H{t),H{t'))(H{t)+H{t'))(H(t) + H{t')-l){t' ~t) H{t)+H{t ' ] ' 2 if H(t) + H (?) + 1 



(7.10) 



Finally, define also for < t, £' < 1, 



A(*,t'):=^^(ff(*),ff(f),*,* / ) (7.H) 



dHdH' 

\2 



_ x (2m-2)!(E?= 



^ ^1^2 l -S-L{H(t),H(t')) sin (/?(#(*')) - /?(#(*))) if ff(t) + H{t') = 1 

B OT (f,i):=S \2 2m-i (7.12) 

(-l) m+1 2 :° ml) 2 L 12 (H(t),H(t')) J] (#(*)+ If (t')-f) ifff(t) + H(f)^l, 

with /? = if a + = a~ . Remark that A(t,t') and B m (t,t') are bounded on (0, l) 2 . 

Lemma 4 Under Assumptions of Proposition^ and for any < e < 1/2, [ne] < k < k' < [(1 — e)n], k' — k > 2q. 

cov{y:x {M ^),v^x {a+ ^)) = (v:h£),v:-) x *^ x ( K V^)x a -) T 

k k' \k' - k\ Hi -^ +H ^- 2m 
+ l {m > 2} B m (-, — ) ^TTTTItt + A„(fc, fc'). (7.13) 

The remainder term \ n (k,k') in satisfies the following bound: for any 5 > there exist no G N*, fco € N* suc/i 

£/ia£ /or any n > no, |fc — fc'| > ko, [ne] < k < k' < [(1 — e)n] 

/ 1 I _ U\H(k/n)+H(k' /n)-2m . 
M*.*0I ^ ' + ' n H( fc /n) +g (i,/n) )• 

Remark 10 In (17.1 3|) V,f — = — (S iaA or depending on whether m = 1 or m > 2 holds. In particular, for a = a* 

n n ^ — ' 

i=l 

the quadratic form on the r.h.s. of f|7.13[) coincides with the term ^-) V^*iJ(^) V„ a *iJ(^-) in the expansion (|4.6[) . 
Proof of Lemma^ For fe, k' € {1, 2, • • • , n — q — 1} and i, j <E {0, 1, • • • , q} denote 

5 i ;=H(- + -)-H(-), 5>:=H(- + l)-H(-), \\6\\ := max \SA, \\S'\\ := max $|. (7.15) 

Assumption -ff(-) € C(0, 1) implies that there exists C > not depending on k and n and such that max{|o~i|, \ d'j\} < 
Cu"' ,A1 . Using the definition in (|7.ip and the notation in (|7.15l) . 

E[v;r^(a + ,a-)(-)v; a A (a+!a - ) (-)] = £ a, aj q(h(-) + s it ff(-) + - + -, - + 3 -) . (7.16) 

i,3=0 

Similarly to (|7.9p . for any collection p — {pi, ■ ■ ■ ,Pa) € N 4 of integers and s^s', introduce partial derivatives of order 
|p| :=pi H hp 4 of Q in (J71J: 



«u-,^3,^ c»y) - dHP1 jr'L^ (g(s) ' g(s,) ' s ' s,) - (7 ' 17) 
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Without loss of generality, assume k' — k > 2qm the sequel. Using the Taylor expansion of order 2m of Q in (H, H' , s, s'), 
one obtains: 

E[^ (0+ja -)(^)v;«jf (a+i0 - ) (^)] = j -(p)+ E ( 7 - 18 ) 



where 



0<|p|<2m |p|=2m+l 



^ ^ i,j'=0 



up) < c sup K?] ) ... iS , P4 ( S , S ')i e ic^p^n-n 

j{ k k' \ . . ' n n 

and where the set i(|, ~) C (0, l) 4 is defined by 

/(*, -) := {{H, H', s, s') : \H < - < |« - -| < «, |*' - -| < « , « < A 

n n I ' n 1 1 n 1 1 n 1 n ' n ' n ) 



Note terms in (|7.18[) corresponding to \p\ — and \p\ = 1 vanish because of the moment condition (|2.1[) . In particular, 

,k k\ j k\ 



E ^(p) =E a * (* (-.-)+-+ ' ^ -)) = 0. 
^— ' V n n n n n ' 

\P\=X i,j=0 



In a similar way, 

T,Mp) = ^2(E^)(E«^e'(^) + --' + 2(E\)(E\)€(" 

|p|=2 i=0 i=0 ' i=0 j=0 

= ( W *) > ^*)xiiP)(^ I ^)x(^(^) > ^) T ) (7.19) 
n n n n n n 

since Yli=o a A = Kf-^(^): ' ' ' >J2j=o a ji = Note that the partial derivatives in (|7.17|) involving arguments U 

and £f' contributes to a factor log ( fc ~ fc ) at each differentiation while each differentiation in s or s' contributes a factor 
I k ~ fc I : in other words, 



sup 



( S , S ') < Cfl+l^l^^^-^xllogr^ir^). (7.20) 



Also note 

\V?H(-)\ < Cn'^ Am \ |V„ a -| < Cn- rn l<k<n-q. (7.21) 
n n 



From (f7T2T1) and (f7T20|) we obtain 

^ 1 fc |g(fc/n)+g(fc7n)-2 . 

| 2^ ^ 6 V ?7 2(r ) Am) H n H(k/n)+H(k> /n) 1 {m=l} J 

|p|=2 

since F„ a ^ = for m > 1. We also have from (f7T2"j) (with /3 = if a+ = a - ), 
J„(0, 0, m, m)l {m > 2} = (-l) m+l (2W " 2) ' ( ^°^ m)2 £(g(-),g(-)) sin (fi{H&) - /?(#(-))) '* " * 



(7.22) 



/ll— 2m 



2(m!) 2 n'n n n n 

ifff(*)+ff(£) = l; 
(T q a i m ) 2 k V 2r ^ 1 k V \k - k'\ H ^)+ H ^)- 2 ' m 



2(m!) 2 v V' 11 v v n y y n' ' „h(±)+h(%-) 

£=0 11 



ifff(£)+tf(£)^l (7.23) 
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which coincides with the second expansion term on the r.h.s. of (|7. 13|) . It remains to prove that all the other terms in 
the Taylor expansion (|7.18p are negligible with respect to the main terms in (|7.13[) . in other words, to show the bound 
(|7.14[) on the remainder term X n (k,k'). Define 

log 2 n \k' - k\H(k/n)+H(k' /n)-2m 
A •= - u J ! 

™ ' n 2(rjAm) ' n H(k/n)+H(k' /n) 

We shall prove below that there exist fix > 0, /3 2 > independent of n, k, k' and such that for any < e < 1/2, n e N* 
and [ne] < k < k' < [(1 - e)n], 

J n (p) < S n A n , 2<|p|<2m, p^ (0,0, m,m), (7.24) 
Jn{p) < S' n A n , \p\=2m + l. (7.25) 

where 



S n (k',k) 



( fx + (J -) f , S' n (k', k) := <5„(fc', fc) + — — . 

n n \k — k\ 



The statement of the lemma including (|7T4|) follows from ([7718]) . (|7T9f and (fT2i|) . (|7T25j) . 
Let us prove (|7.24p . Define 



ET B (p) := 5:^(1)", 0» := $>(3> Pl (£r- (7.26) 

i=0 3=0 

We claim that for any 77 > and m > 1 

|fn(p)| < Cf^— + ^^) (7.27) 

and a similar bound holds for f/^ (p) . Let us check that ()7.27|) and (|T. 20[) imply (|7.24|) . Indeed, from the above relations 
and the definition of J n (p) we obtain 

\Jn( P )\ < |^K..., t ',4(^)|l^(p)ll^(p)l 

< cao&n^+p- i fc, ~ fc i g(fc/ " )+g(fc,/ " ) - p3 - p V 1 I — V — I 1 

— \ £> J I n I \ n mVp 3 n V+P3 J \ n mVp 4 n V+Pi 

< C (lognf 1+ P 2 (7 U + 7 12 + J 21 + J 22 ), 



where 



_ <k' — k iH(k/n)+H(k' /n)—p3—p4 1 1 fc' — iH(fc/n)+H(fc'/n)-ps-P4 1 

7 11 : ~ 9^7: ^12 ' 



n ' n 2m ' I n I n (mVp4)+i)+P3 ' 

_ I fc' — |H(*/n)+H(fc'/T»)-ps-p4 1 _ \k' — k \H(k/n)+H(k' /n)-p3-pi 1 

-'21 •— ~ T^JZXZ^J^T. i ^22 •— 



n (mVp 3 )+?7+p4 ' 1 rt 1 n 2»7+p3+p4 ' 

It suffices to prove (|7.24[) for 7j, i, j = 1,2. 

Let H := (H(k/n) + H(k'/n))/2, j = \k' - k\, p := p 3 + p 4 . For 7 U , (gill follows from 

(-) ^r = y (-)^r^(-) (7) ^ = (^jH(-) < a„<W *n,»:=-. (7 ' 28) 

which is immediate from < p < 2m — 1 and 1 < j < n. 

Next, consider 7 2 2- Assume first that rj < m. Observe that in this case, A n > rr 2r > if j > n('? _ - H V( TO_ - E and 
A„ > j 2H -^ m n - 2H if j < n (n-H)/(m-H) _ Thug in the cage ^ > n ( v -H)/(m-H) t the bound (fT24j) for 7 22 translates to 

(-) SUI - p m < with W-T = 4(-) 2H <(-) 2H ( 7 - 29 ) 

n n zr i+p n zr i n ztl J p n n 



since j > 1. Next, let 1 < j < n^/M.p > 1, then the bound (f77M|) for 7 22 follows by 

j2H-2m j2m~p 

„>„+,. "^77-- with fe ." := ^r 



(^) 2 ^^ < ^^T^' with = (^r) 2m ^- fe <n-fe, (7.30) 



26 



where P22 '■= 2H ^ m JlMMs. H *> > 0. Next, consider I22 for rj > m. Then the first term in the definition of A„ is negligible 
with respect to the second term and the corresponding relation reduces to 

■2H-p -2H-2m -2m-p ■ 

< S22, J 2ff , With fa tn ;= J —— < (^) 2m . (7.31) 



n 2 V +2H - n 2H ' ' n 2 ^ 

Next, consider 1 \ 2 . Accordingly, we need to show 



A ryZJ 1 1 A.2H ~2m 

( J_\ 2H -P3-Pi 1 < <: ( 1 I J \ (7 19) 

\ n > n (mV P4 )+>)+P3 - ™V n 2(^Am) + n 2H )> K 1 -^) 

where S n satisfies the bound in (|7.26[) . Let r\ < m then (|7.32l) becomes 

A „ „ 1 1 -2H-2m 

\ n ) n(mVPi)+ri+P3 > ° n \ n 2 V + n 2H >■ y-^) 

Let j > n™- H . The first term on the r.h.s. of (|7.33l) dominates the second one and (|7.32[) for 2H < p 3 + p 4 follows from 
,U2H- P3 - Pi 1 < s , 1 with S[ 2n :=n-^, (3' 12 >0, 



where 



# 2 := 2H-T 1 +(mV Pi )-p 4 + ^—^-( P3 + Pi -2H) 

m — H 



(m- H)- 1 



(.P3 + P4 ~ to)(t/ — H) + H(m — rj), p 4 > m, 
(to — pi + H)(m — rf) + p 3 (rj - H), p 4 < to. 



Thus P' 12 > H(m - rj)/(m - H) > 0. On the other hand, if p 3 + p 4 < 2H < 2 then the first inequality in (|7.33[) holds 
with S n = 5' 12 n where 

X> — (3\ 2H -P3-P4 1 < „-(m-r,) 

with to — rj > 0, since m > 2, p 4 < 2 and (m V p 4 ) — r] + p 3 > rn — r\ > 0. 

Next, let j < n™- H . Then the second term on the r.h.s. of (17.33)) dominates the first one and (|7.32[) follows from 

A „„ i a2H— 2m „-2m— p 3 — p4 

/J\2iT-p 3 -P4 1 . c// J -,i if// _ J -/3" 

n (mV P4 ) +I7 +P3 S ° 12 >" n 2H ' WltJl °12,n — n „+(mVp 4 )-P4 S " 

where 

// r7 — if 

/3" 2 := J) + (roV])4)-p4-( -)(2m - p 3 - p 4 ) 

m ~ H 



(m- H)- 1 



(P3 +P4 — to) (7/ — + _ff(m — 77), p 4 > m, 
(to — jj 4 + i?)(m — rf) + p 3 (rj — H), pi < to., 



and therefore f}'{ 2 = f3'i2> implying j3" 2 > as above. 

It remains to consider I12 for 77 > to. In this case the second term on the r.h.s. of (|7.33[) dominates the first one for 
ah 1 < j < n and (|7T32| follows from 



A „„ 1 -2H-2m A2m-p 3 -p 4 ■ 
'J\ 2H -Ps-Pi 1 < c/// J s-m ._ J_ K (l\2m-p 3 -P4 < (I) 



since p 3 + pi < 2m and rj + (to V p 4 ) — 2to + ^3 = (77 — to) + ((to V p 4 ) — to) + j»3 > 0. This proves (|7.32p or (|7.24|l for 
Ji2- Since consideration of I21 is completely analogous, the proof of (|7.24j) is now complete. 

Let us prove the claim (|7.27l) . Write 5i = 5 0i + Ti, where for i £ {0, • • • , q} and r = {0, • • ■ , [77]}, 

W IT ' t. 
= ^f^ (1) r Hr : = H(r} ( *). 
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By Assumption 0(77), 
From binomial expansion, 

U n (p) 



max(|<5 0i |) < Cn' 11 , max(\H r \) < C, 



r H H"[»,]=Pl 



Pi 

T 0, ■ ■ ■ >T[ V ] 



[v] H T -\ q i 



r=l 



i=0 



where t% € N, i = 0, 1 • • • , [rj\ are nonnegative integers. According to (|2.ip . for any s = 0,1, 

. 9 



i=0 



Va,(-) S < Cn- S l(s>m). 



Therefore from (|7.34[) it follows that 

i=0 

Note that for tq = 0, we have J2r=i rT r ^ Pi- Hence 



i7-o/l-\P3+Er=l r " r r 



0, 



r = and p 3 + ^[=1 rr r ^ 



otherwise. 



lv] 



P3-Elii^i( T = o,p 3 + 2 r ' T '->«i) < (n _m )A(?i- pi - p3 ) < 



-(mVpa) 



Therefore 



< C(n -,?-P3 + n~ (mVp3) ), 



i=0 



(7.34) 



proving (|7T27|) and (|7T24|) . 

It remains to prove (j7^5|) . Let L/„(p) := £? =0 Mf(£) P3 |- Since ll^ll ^ C^ 7 ^, see (17341 . it follows that 

U„(p) < Cn-^)n-p 3 (7.35) 
and a similar bound holds for U'Jp) := EL | a i( 5 j') P2 (£) P4 |- Relations (TT351) and (TT201) imply 



J»(p) < qi 



k' - k,2H-p 



x I log (*L h\ |2™+1-P\ n -(r,Al)(2m+l-p)-p 



(7.36) 



where p = p?, + pa £ {0, 1, • • • , 2m + 1} and \p\ — 2m + 1. Relation (j7.25[) now can be verified directly, by inspecting 
the cases 77 < 1, rj > l,p = 0, 1, 2m + 1 and 1 < p < 2m + 1 separately. This ends the proof of Lemma [4] □ 



Proof of Proposition QJ Note first that the result of Lemma Q] cannot be used for establishing Proposition Q] because 
Proposition Q] treats the case k close to k! which is excluded by the proof of Lemma 0] (the second derivative R t ,t' does 
not exists for t = t'). 

For any [tie] < k < [(1 — e)n] similarly to (|7.16l) and using the notation in (|7.15[) with k' = k + 1 we rewrite the variance 
V n (k) :=Var(y n a X (Q+ , a - ) (^)) as 



Vn(fc) = \ V a i a j Q(H(-) + 6 i ,H(-)+8 j ,- + -,- + t). (7.37) 
2 ' \ n n n n n n/ 

i,j=0 



Next, we use a Taylor expansion of Q given by two analytic expressions in (|7.2j) . see also (|7.8[) . In order to avoid 
cumbersome notation, we shall assume in the rest of the proof that the function H(-) is separated from 1/2, i.e., 
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sup te(01) \H(t) - 1/2| > 0. Then the second expression for Q in (J772J) applies, viz., Q(H,H',t,t') = (±)(Q n (H, H',t) + 
Q 22 {H, H', t') - Q 12 (H, H', t, t')), where 

Q n (H,H',t) :=L 11 (H ) H')t H+H ', Q 22 (H,H',t) := L 22 (H, H')t' H+H ' , Q 12 (H, H> ,t,t') := L 12 (H,H')\t' - t\ H+H ' . 

Accordingly, V n (k) = |(V n ii(fe) + V n22 (k) — V n i 2 (k)), where V n %j are defined as in (|T.37[) with Q replaced by Qij, 1 < 
i < j < 2. Note that (H, H' ,t) £ (0, l) 3 M- Qu(H, H' , t) is analytic except for H + H' = 1 and the last case is excluded 
by our assumption that H(-) is separated from 1/2. Using the last fact, the Taylor expansion of Qn as in (|7.18l) in the 
proof of Lemma 2] and the bounds (|7.21[) . it follows that for all n large enough, 

sup \ Vnll {k)\<C{e)^— (7.38) 

with C(e) > not depending on n. Clearly, the same bound (|7.38p holds for V n22 . 
It remains to estimate V n \ 2 . For i, j G {0, • • • ,q}, we have 



expUdi + d,)log f = +0 



Here and in the rest of the proof, we denote O(-) a bounded function depending on H (•), e and q only, and not depending 
on k or n. Next, since (H,H') i— > Li 2 (H, H') is analytic for i7 + _ff ' 7^ 1 , so 

Wtf(-) + * , ff(-) + = 1 + o(^ T ). 

Therefore, 

V„ 12 (fc) = g(g(fe/n), a) n~ 2H ^ + J ( °| } ^ A1 ) (7.39) 

where C(H,a) = —n 2H Y^ij=Q a i a j |^^| 2 = 2Var (V a P>if (0)) does not depend on n and does not vanish for H E 

(0,1), see Sec. 2.2. 

Using (jTUgl) and (TH?9")) . we obtain: 

Vn(fc) = --C(g(fc/»),a)n- 2H ^)+of J A +0( ^ ) 

= Var(U n a B„ (M (-))(l + 0(^)+o(— J —)). 

Hence the bound of Proposition [T] can be deduced. □ 



Proof of Corollary Fact (i) is immediate from Proposition Q] By Property [1] E(V^ X( a + a -){k/n)) 2 ~ (4 — 
A H(k/n)^ n -2H(k/n) and hcnce E(U„ Q * X (a + , a - } (fc/n) ) 2 > Cn- 2H{ - k ''^ for some C* > independent of 1 < k < n - 2 
and n. Use the last fact, Proposition [2] gU), gj}, boundedness of A(t,t') and B(t,t'), and |U,f < CrT^ 2 ) to 

conclude inequality (|4.8I) for \k — k'\ > k and fco large enough; for \k — k'\ < ko (14.81) is trivial since the l.h.s. of (|4.8I) 
is bounded by 1. 

It remains to prove fact (ii). Let k' — k =: k > 1, H ;= H(k/n), H' := H(k'/n). Let first k < n 1 ^". Then in view of 
(|4.8I) it suffices to check that 

The second inequality in ([7T4H|) holds for k < ti 1 "" because of JOJ and min(2(r7 A 2) - 2F, 2(7? A 2) - 2i?') < 2(77 A 2) - 
H - H'. The first inequality in (|7T40"| holds trivially since 7 - (4 - H - H') < 0. Next, let > n 1 "". Then using (|4~8f 
it suffices to check the inequalities 

1 C 1 C 

yA-H-H' - n (7-0)(l-a)' n 2( V A2)-H-H' - n ( 7 -S)(l-a) " 
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which easily follow from (|4.9p . similarly as above. 



□ 



Proof of Proposition [3J . Part (i) is immediate from Theorem [T] (i) and Corollary [3] (i) . Part (ii) follows from Theorem 
Q] (ii) and Corollary [3l by taking 9 — and 7 > 1/2 sufficiently close to 1/2 (in such case the inequality 7 < 4 — 2H(t) 
of (|4.9j) is satisfied since H(t) < 1). To prove (hi), use Corollary [3] and Theorem Q] (hi) with = and 7 > 1/2 
which can be chosen under (|4.10j) . with [i = (2 — 3a) /4 and condition (|4.11|) . Note that the lower bound (|4.11|) implies 
(77 A 1) > (2-3a)/4 and (jlTO)) implies 2 ((77 A 2) - sup te(01) H(t)) > (2-3a)/4, or condition k > /i of Theorem Q] (hi). 
Parts (iv) and (v) follows similarly from Corollary^ Theorem [1] (iv) (with 9 = 0, k > [i\ = (1 — 3a)/4 satisfied from 
(|iT0]) ) and Theorem [2] (with a > 1 - 2k and (1 + 2(77 A 2))" 1 > 1 - 2(77 A 1) for all 77 > 0). □ 
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